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ABSTRACT 

Safety analyses of a loss-of-coolant accident in a 
pressurized-water-reactor nuclear power plant require 
modeling the flow through the main-coolant pumps in any one 
of several extreme operating conditions. A semi-empirical 
method of predicting the two-phase performance of centri¬ 
fugal pumps has been proposed and was applied to normal pump 
operation with good results. In this report the proposed 
method is applied to the reverse-flow, reverse-rotation mode 
of operation, also with good results. 

By making simplifying assumptions, the flow through 
the pump is idealized and expressions for the theoretical 
pump characteristics derived from Euler's equation for the 
mean flow through a turbomachine. General expressions for 
the flow losses in two-phase flow are derived, and the flow 
losses expressed as the difference between the theoretical 
characteristics and the actual performance of the pump. 

The ratio of the losses in two-phase flow to those in single¬ 
phase flow is shown to be primarily a function of pump 
geometry and fluid void-fraction. Using experimental data 
from single- and two-phase flow tests on a pump of known 
geometry, the ratio of the losses was calculated for various 
conditions of flow and plotted as a function of void fraction. 
The results plotted in this fashion show excellent correlation 
and the function can be used to predict the two-phase 
characteristics of similar pumps. 
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LIST OF SYMBOLS 


A 

flow area 

a 

P 

Ot ^ V 

two-phase flow function = [t-] — 

l-a 0^ 

c 

fluid velocity, absolute 

^tp 

two-phase flow function (eq. 8a) 

g 

gravitational acceleration 


constant in Newton's Law 

H 

head 

H* 

AH ^ - AH ^ 

head-loss ratio s t- ^ ^° ^ 

''“ospth ■ 

h 

enthalpy 

K, K^, 

K^, Kg constants defined by eqs. 12, 13, and 15 

• 

in 

mass flow per unit time 

N 

rotor speed, rev/min 

P 

fluid pressure 

Q 

volume flow per unit time 

r 

impeller radius 

R 

1-x 2 

friction-factor multiplier = 

s 

slip velocity ratio = (C /C^) 

U 

rotor peripheral velocity 

w 

fluid velocity relative to rotor 

X 

quality = (m^/m^) 

a 

void fraction = (A^/A^) 

B 

angle of fluid vector relative to tangent to 
rotor periphery 

B' 

angle of tangent to rotor-blade centerline relative 
to tangent to rotor periphery 
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property difference 
deviation angle 
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e flow ratio e ^^be^^tp^ 

p fluid density 

(b flow coefficient E (C /U) 

^ m 

2 

ill work coefficient = (g Ah /U ) 

CO 

2 

ip' head coefficient = (gAH^/U ) 

SUBSCRIPTS 

o total (static plus dynamic) or stagnation value 

of property 

1 plane at entrance to rotor 

2 plane at exit from rotor 

3 plane at entrance to pump 

be best-efficiency point 

f wall-friction component 

L liquid 

m meridional or radial component of velocity 

6 tangential component of velocity 

s flow-separation component 

sp single-phase component 

T total flow, liquid plus vapor 

tp two-phase flow 

th theoretical or ideal conditions 

V vapor 
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INTRODUCTION 


BACKGROUND 

In a nuclear power plant of the pressurized-water- 
reactor type (PWR), sub-cooled water at high pressure and 
temperature is circulated through the reactor core by several 
main-coolant pumps, each in its own coolant leg. Safety 
analyses of postulated loss of coolant accidents (LOCA) 
require models of the flow through the main-coolant pumps 
for any one of several operating modes in both single- and 
two-phase flow. In an effort to more precisely predict the 
contribution of the main-coolant pumps to flow through the 
reactor core during a LOCA, there has been a continuing 
effort to improve these models, particularly for two-phase 
flow. 

A semi-empirical method for predicting the two-phase 
flow performance of a centrifugal pump was proposed by 
Mikielewicz and Wilson,^ and was applied to normal pump 
operation (forward rotation and flow) with promising results. 
Details of their procedure are summarized below. 

Starting with the Euler equation a theoretically- 
ideal head-versus-flow relationship is derived for 
both single- and two-phase flow. 

By examining the hydraulic losses in two-phase 
flow it is shown that for a given pump the relative 
pressure losses—the difference between theoretical 
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and actual pressure rise—are related to the relative 
pressure losses in single-phase flow by a function of 
void fraction, a property of the two-phase fluid. 

The ratio of relative pressure losses in two- 
phase flow to those in single-phase flow is defined 
as the "Head-loss ratio". 

The head-loss ratio is then calculated by 
combining the derived theoretical characteristics 
with experimental data for single- and two-phase flow. 

For first quadrant operation, good correlation is 
shown between the head-loss ratio and void fraction, 
and the result can be used to predict the two-phase 
characteristics of the pump. 

It is shown that the head-loss ratio for a given 
pump is related to other pumps by a function of pump 
geometry; further experimental work is required to 
determine this relationship. 

PURPOSE 

The purpose of this report is to determine if the 
procedure is applicable to the third quadrant of operation 
(reverse-flow and rotation). The steps of the procedure, as 
outlined, will be modified only as required by the different 
flow geometry, and the terminology used will be the same. 
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PROCEDURE 


DERIVATION OF THEORETICAL CHARACTERISTICS 

A centrifugal piomp subjected to conditions such that its 
normal flow and rotation directions are reversed will operate 
in a manner similar to a radial-in-flow turbine. The Euler 
equation can be applied to determine the characteristics of 
the pump in this mode of operation: 

gc = Uj^Cg^ - ^2^Q2 

where: 

Cg = fluid tangential velocity 
g^ E constant in Newton's Law 
Ah^= change in total enthalpy = hol''^o2 
U = rotor peripheral velocity 
subscript 1 E rotor entrance 
subscript 2 = rotor outlet 

The fluid tangential velocity at rotor outlet, 
could be significant because conditions of flow and rotor 
speed would be other than optimum and because, in general, 
a pump is not designed to be an efficient turbine. For 
these reasons the simplifying assumption of zero outlet 
swirl will not be made. This general form of the equation 
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is exact and will yield the actual enthalpy change across 
the pump, for any flow condition, if the proper values for 
the velocities are used. 

In the inlet region of the pump. Figure 1, precise 
flow angles and velocities are unknown and an assumption 
must be made relating to the flow rate entering the 

scroll. Using this assumed value of Cg^ the velocity tri¬ 
angle at the rotor entrance, as shown in Figure 2, can be 
solved for all velocities in terms of the flow rate. 

At the rotor outlet, flow geometry is specified by the 
design of what is normally the impeller inlet eye. A typical 
design with swept-back blading is shown in Figure 3. Reverse 
flow and rotation will produce a velocity diagram as shown in 
Figure 4, which can also be solved for all velocities in 
terms of the flow rate. 

If the reverse flow through the pump is idealized by 
making simplifying assumptions, the Euler equation can be used 
to approximate the theoretical head and flow characteristics 
of the pump. The following flow conditions are assumed. 

1. That the fluid tangential velocity at rotor entrance, 
Cgj^, is equal in magnitude to the absolute fluid 
velocity, C^, entering the scroll. 

2. That the relative flow angle at rotor outlet, B> 2 ' 
is constant for all flow rates and conditions. 

That the deviation angle at rotor outlet, 62 » is 
zero, that is, 62 =^ 2 • 
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PLANE 3 



FIGURE 1 - TYPICAL CENTRIFUGAL PUMP IN REVERSE-FLOW AND 
ROTATION 
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FIGURE 2 



Uj^ - rotor peripheral velocity 

- fluid absolute velocity 
Cg^- fluid tangential velocity 

fluid radial velocity 

- fluid relative velocity 

- fluid relative flow angle 

- VELOCITY DIAGRAM AT IMPELLER INLET 
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FIGURE 3 - TYPICAL IMPELLER-EYE WITH REVERSE-FLOW AND 
ROTATION 
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FIGURE 4 
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- VELOCITY DIAGRAM AT IMPELLER OUTLET 
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with the above assumptions, equation 1 becomes 


Ah 


o th 


= 


- U2S2 


( 2 ) 


where 

subscript th e theoretical or ideal conditions 

✓ 

For the specific case of loss-less incompressible flow, 
the theoretical total enthalpy change is proportional to the 
change in theoretical total head.^ 

change in total head = 
gravitational acceleration 

and the expression for the theoretical total head character¬ 
istics of the pump is: 

5 “cth = “iSl - “2=92 


where 

AH^ E 

g = 
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Theoretical single-phase-flow characteristics 


For single-phase flow, the velocity triangles are solved 
by using the continuity equation and the assumptions made to 
idealize the flow. 

In the pump inlet region, it was assumed that 


C 


C 


ei 


3 


where 


subscript 3 = plane at inlet to pump 

From continuity 


where 

Q = volume flow per unit time 
A = area normal to the flow 

subscript m = meridional or radial component of velocity 

therefore, 



( 4 ) 
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At the rotor outlet, from Figure 4, 


'62 


= U - _ 

^2 tan 


Since 


U 


2 



r 


2 


r 


1 


where 


2 ... 
— = rotor-outlet-to-inlet radius-ratio 


and by continuity 



A 

^ml A 


1 

2 


then. 


C 


62 




ml 

tan B, 


(5) 


By substituting equations 4 and 5, equation 3 becomes 


^ ^^othsp ■ ^l^ml ^ 



fna__, 

tan $2 


where 


subscript sp = single-phase flow 
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Dividing by and defining the head coefficient, as 


gAH 

'i'i = -T 

uf 


and the flow coefficient, as 


_ 


U. 


equation 3 becomes 


^9 9 ^ 9 

♦ithsp = -<?7> ^ 'a 7 ^ J- tan 6, ’ 


( 6 ) 


This expression is a linear function of ip ' versus (p , 
and is similar to the expression derived for first-quadrant 
operation. 


Theoretical two-phase-flow characteristics 

In two-phase separated flow it is assumed that the 
liquid and the vapor phases have separate mass flows and 
that their velocity components may also be different. In 
this case, equation 2 becomes 


U. 


■^ithtp ^ ^^“^i^Sli'^^iS vl^ 


U. 
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where 


• • 

X = quality = 

and m = mass flow per unit time 
subscripts L = liquid 

T = total flow, liquid + vapor 
tp= two-phase flow 
V E vapor 


From the velocity triangles at rotor inlet and outlet 
and from the assumption that Cq^ = C^, 


’^Ithtp 


1 ^1 ^1 
— [(l-x-)C TT 7— + x-C , 

1 mLl 1 mvl 


U" C . „ 

^ r / n \ / 1 mXj ^ 

^[(1-X2) (1~U2 tan 


)+x (1—_)1 

’ 2 ^ tan ^ 


From continuity 




L 


(1-x) 

Pl \ ” Pl(i-«)a^ 


m 


Snv 


V 


X 


P A 

V V 


0 a A 

T 


where 

a = void fraction 
p E fluid density 
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and by defining the two-phase-flow coefficient, 




Qj. rn_ 

= _tE = JL 


tp- A U Ptp ^ ^ 


where 


= mean density of two-phase flow 

= (l-a)Pj^ 


Equation 2 becomes 




2 PlI , 

X, - - + X 


‘1 a. 


vl 


2) - 

1 uf ^2 


a. 


1 +- 


r, 2, 2^ '^L2 2, ^ 

[(1-X-) + — - -- X ]} 

2 a p 2 


_ ^v2 . 

(l-'^2) ^L2 


(7) 


By defining a two-phase-flow function, a. 


a = 


a 


^1-a^ 



and a two-phase slip-velocity ratio, s. 
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The quality of the fluid can be expressed as 


as 


By substituting these expressions, equation 7 becomes 


; . - * (1+aj^) (1+aj^Sj^) ^ r2 2 

Ithtp tpl (i+a^s^)^ ^1 


[ 1 - 


<P 


tp2 


(l+a^) (l+a^s^) 


tan So /1 \ 2 

2 (l+a2S2) 


-] 


( 8 ) 


Applying the continuity equation to the two-phase flow 
coefficients at rotor inlet and outlet yields. 


p. T A, U. 

A - A tpl 1 1 

♦tp2 - ♦tpl 


which reduces equation 8 to 




♦tpl ^1 "ij 


Ithtp A^ tpl tpl tan 62 tp 2 P^p2 ^2 ^2 


where 


tp 


= (1+^)(1+as ) 
(1+as)^ 


(8a) 
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Rearranging: 


w . - 2,2 , r 1 . , -^2 ^tpl "1 _tp2_, , 

^Ithtp rj_ tpl P^p2 A^ tan 62 '^tpl 

( 9 ) 

If the flow becomes single-phase, equation 9 reduces to 
the single-phase case, equation 6, because ^tp2 

P^p^/Ppp2 all equal unity when the void fraction, a, 
equals zero. 

As derived, equation 9 is cumbersome and requires a 
knowledge of the relationships between inlet and outlet 
conditions to calculate values of the head coefficient for 
various conditions of void fraction and flow. By assuming 
that fluid properties at rotor inlet and outlet are equal, 
equation 9 reduces to 


'^Ithtp 





^1 ^2 ^1 

^aJ A2 tan 82^ ^tpl '^tpl 


( 10 ) 


The assumption is valid in this case because the moment of 
momentum at rotor outlet, U C _, is very small in comparison 
to U.C-T, and the outlet properties appear in the term derived 

X 91 

from ^20^2' thus the effect of changes in fluid properties 
from inlet to outlet is minimal. 
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The two-phase flow function, f 


was defined by 


tp' 


equation 8a as 


tp 


(1 +a)(1 +as ) 
(l+as)^ 


where 


a s 


(—) — 
^1-a^ Pl 


s = 


V 


For a steam-water mixture, liquid and vapor densities 

are given by steam tables for various pressures. The slip- 

velocity ratio, s, for two-phase flow in a horizontal pipe, 

2 

was correlated as a function of pressure and void fraction 

as shown in Figure 5 . A correlation more applicable to two- 

phase flow in centrifugal pumps was not available; therefore, 

the relationship described, though for vastly different 

conditions of flow, was applied. 

Values of the two-phase flow function, can be 

calculated as a function of void fraction, a, for various 

pressures, as shown in Figure 6. 

The theoretical characteristics for third-quadrant 

operation of a centrifugal pump in both single- and two- 

phase flow can be represented as shown in Figure 7 . In 

single-phase flow, f^^ is unity; in two-phase flow, f^^ is 

determined by pressure and void fraction. 
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FIGURE 5 - SLIP-VELOCITY RATIO, s, VERSUS PRESSURE 


(from reference 2) 
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FIGURE 6 - TWO-PHASE FLOW FUNCTION, , VERSUS VOID 

tp' 

FRACTION, a, FOR 200 AND 500 psia 
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FIGURE 7 - THEORETICAL HEAD-VERSUS-FLOW CHARACTERISTICS 
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FLOW LOSSES 


At a given rate of flow through a centrifugal pump, 

the difference between the measured head and the theoretically 

ideal head is the result of flow losses generated within the 

pump, and changes in the exit deviation angle, the 

assumed ideal direction. For third-quadrant operation, in 

single-phase flow, this difference is shown in Figure 8. 

The curve of actual head versus flow is typical for this 
3 4 

quadrant, ' and the theoretical relationship is approxi¬ 
mated from the previously derived expression (equation 6). 

The actual head, AH , is a function of the flow-rate, Q, 

o 

squared;^ therefore, with flow increasing from zero, the 
losses decrease as the square of the flow^ pass through a 
minimum, then increase as the square of the flow. The flow- 
rate corresponding to the point of minimum losses is de¬ 
fined as the best efficiency flow, third-quadrant 

operation. 

In the first-quadrant of operation, the flow losses 
a 

behave in^similar manner and the expressions derived to 
relate the losses in two-phase flow to the losses in single¬ 
phase flow may be applied to the third-quadrant. The 
constants and flow-rate ratio are those applicable to the 
third quadrant. 


28 




FIGURE 8 - TYPICAL THEORETICAL AND ACTUAL PUMP 
CHARACTERISTICS 
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The flow losses are expressed as the sum of wall-friction 


losses and flow-separation losses, 

Ap = Ap - + Ap 
^o “^of “^os 


( 11 ) 


where 

Ap^ = loss of total, or stagnation, pressure 
subscripts f = friction component 

s H flow-separation component 


The wall-friction component of the losses in two-phase 
flow is given by 


Ap , ^ = R K 

•^otpf 


;^)2 


Q 


tp 


( 12 ) 


where 

5 1-x 2 

R = friction factor multiplier 

K = constant at high Reynolds numbers; includes pump 
geometry, the single-phase friction factor, and 
liquid density 

By approximating the flow-separation losses as proportional 
to the square of the difference between the existing flow 

g 

and the flow at best efficiency point, the two-phase 

flow-separation losses become 
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(13) 


Ap 


otps 


= K. 


•tp 


[(l-a) ( 


1+a 

1+sa 


) - e] 


where 


= constant, depends on flow regime 
" - °be/°tp 


Adding these components yields the expression for the 
total losses in two-phase flow, Ap^^^, 


''Potp = ^ 'Ff'' °tp * "4 °?p'<l-“> 


substituting 


^tp _ rl“Ce, rl+a , 

4-x^ ^1+sa^ 

The total losses in two-phase flow become 


Ap 


otp 






(15) 


where 

and Kg are constants which include geometric and 
flow-regime factors 

for single-phase flow, equation 15 reduces to 
'^Posp = Kgd-O^l q2 
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The comparison between single-phase and two-phase flow 
is made at like values of flow-rate; therefore, the ratio 
of two-phase to single-phase losses is 


Ap 


/1+a \2 „ r/i \/lta . 1 ^ 


osp 


1 + (l-£) 

D 


(17) 


HEAD-LOSS RATIO 

The total losses, Ap , can be expressed as the difference 

o 

between the actual and the theoretical heads at a given flow- 
rate , 



(AH 

o 


o th 



(18) 


In terms of head, equation 17 becomes 


AH 


AH 


o tp 


o sp 


AH 


o thtp 


AH 


o thsp 



( 1 ±£_) 

'l+sa^ 


1+Kg (l-£)^ 


(19) 


substituting 


tp 


_ ,1-x^ ,1+a 

V T ^ ) V ” 


'1-a 1+sa' 


and X 


sa 

1+sa 
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Equation 19 becomes 


AH 


r 1 + a 


1+a 


otp 


- ^«othtp _ <T+t¥) ^ 


AH^ - AH 

osp othsp 


(l-a) (l+a) [l+Kg(l-e)^] 


( 20 ) 


For the first quadrant, the head-loss ratio, H*, was 
defined as 


H* 


AH 


othtp 


AH 


othsp 


AH 


otp 

osp 


To prevent the introduction of an unnecessary change, 
the first quadrant definition will be used in the third 
quadrant as it is numerically equal to the ratio in 
equation 20. Thus, for third quadrant operation 




( 21 ) 


othsp 


osp 


In this functional relationship: (l+a/l+sa) is a 

function of void fraction and pressure; Kg includes the 
geometry of a given pump and the flow regime of the two- 

phase fluid (a function of void fraction); and e is the 

; 

ratio of best efficiency flow in third-quadrant operation 
to the existing flow. For a given pump, operating at a 
particular pressure, the ratio of two-phase to single-phase 
losses at like values of e is therefore principally a 
function of the void fraction of the two-phase fluid. 
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Substituting the definition for the head coefficient 






gAH 


U 


the head-loss ratio can be expressed as 


H* = ’^thtp ~ ’^tp 

4 ^', - 4 ;' 

^thsp ^sp 


( 22 ) 


CALCULATION OF THE HEAD-LOSS RATIO 

The head-loss ratio as a function of void fraction for 
third-quadrant operation of a particular centrifugal pump 
can be constructed using the derived theoretical character¬ 
istics and the actual characteristics obtained from experi¬ 
mental data. The relationship between H* and a for this 
pump can then be used to predict the third-quadrant, two- 
phase performance of other pumps as single-phase character¬ 
istics are usually known or can be easily determined. 

D.J. Olson conducted experiments on a semi-scale pump 
of known geometry,^ in single- and two-phase flow, and in 
the operating quadrants of interest. Using the results of 
his experiments in the third quadrant, the actual character¬ 
istics of the pump can be calculated. Information obtained 

g 

from the pump manufacturer provided the required constants 
for calculation of the theoretical characteristics for the 
pump using equation 10. 
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The functional relationship between head-loss ratio 
and void fraction for this particular pump is developed as 
described in the following sections. Details of Olson's 
apparatus and test conditions are shown in Appendix D. 

Theoretical head-coefficients 

Details of the semi-scale pump used by Olson are shown 
in Figures 9 and 10. With proper values substituted for the 
constants, equation 10 becomes 

= -0.1896 + 9.155 f^p3_ 

As previously described, this equation may be used to 
determine the theoretical characteristics in both single- 
and two-phase flow. 

Actual single-phase head-coefficient 

A non-dimensional head-flow curve for the semi-scale 
pump can be constructed, as shown in Figure 11, using the 
single-phase data. Curve-fitting the points by regression, 
using the method of least squares, this curve is accurately 
approximated by the expression 

=101.28 + 0.4 
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FIGURE 9 



- DETAILS OF THE SEMI-SCALE PUMP 
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FIGURE 10 - DETAILS OF THE SEMI-SCALE PUMP IMPELLER EYE 
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FIGURE 11 - HEAD-VERSUS-FLOW CHARACTERISTICS OF SEMI-SCALE 
PUMP IN SINGLE-PHASE, THIRD-QUADRANT OPERATION 
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Head-loss ratio 


The head-loss ratio is calculated for each test 
conducted in two-phase flow, third-quadrant operation. 

RESULTS 

The results of these calculations are shown in Table 1, 
and the relationship between head-loss ratio and void 
fraction for the semi-scale pump is shown in Figure 12. 
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FIGURE 12 - HEAD-LOSS RATIO, H*, VERSUS INLET VOID FRACTION, 
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TABLE 1 


CALCULATION OF HEAD-LOSS RATIO, H* 


INLET 

VOID 

FRACTION 

^tpl 

"■tpi 

'^Ithtp 

^itp 

'^ithsp 

'^Isp 

H* 

“l 








0.461 

1.022 

0.030 

0.095 

0.919 

0.089 

0.545 

1. 

0.464 

1.015 

0.018 

-0.022 

0.876 

-0.025 

0.454 

1. 

0.466 

1.022 

0.010 

-0.094 

0.816 

-0.096 

0.419 

1. 

0 

1.0 

0.009 

-0.102 

0.469 

-0.102 

0.417 

1. 

0 

1.0 

0.010 

-0.101 

0.481 

-0.101 

0.417 

1. 

0 

1.0 

0.024 

0.027 

0.570 

0.027 

0.491 

1. 

0.241 

1.015 

0.044 

0.215 

0.824 

0.209 

0.683 

1. 

0.116 

1.009 

0.047 

0.247 

0.811 

0.244 

0.732 

1. 

0.245 

1.007 

0.044 

0.214 

0.848 

0.211 

0.686 

1. 

0.056 

1.004 

0.025 

0.043 

0.629 

0.042 

0.502 

1. 

0.512 

1.067 

0.037 

0.168 

1.004 

0.146 

0.605 

1. 

0.537 

1.023 

0.022 

0.012 

0.956 

0.008 

0.476 

2. 

0.631 

1.033 

0.015 

-0.048 

1.184 

-0.052 

0.438 

2. 

0.401 

1.014 

0.016 

-0.043 

0.788 

-0.046 

0.442 

1. 


81 

88 

77 

10 

12 

17 

28 

16 

33 

27 

82 

02 

51 

70 
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DISCUSSION 


Figure 12 shows that within the range of variables of 
Olson's experiments, the head-loss ratio of the semi-scale 
pump in third-quadrant operation correlates very well as a 
function of inlet void fraction. An average curve drawn 
through the points can be used to predict the two-phase 
performance at other values of void fraction, flow, and 
pressure, up to the maximum void fraction tested. 

At low values of void fraction (up to 20%) the flow is 
largely homogeneous and the losses in two-phase flow are 
approximately equal to those in single-phase. As void 
fraction increases above 20%, the flow becomes separated 
and the relative losses increase sharply, until at the 
maximum void fraction tested (63%), they are more than 
twice the single-phase losses. By the definition of the 
head-loss ratio, its value is unity at a void fraction of 
100%; the characteristics between 63%rand 100% void-fraction 
are unJcnown, but it is reasonable to assume that the head- 
loss ratio will continue to increase then drop sharply to 
unity at the very high values of void fraction. 

The effect of pressure on the correlation of the head- 
loss ratio is not apparent for this limited amount of data. 
Two series of two-phase tests were conducted, one at 200 psia 
system pressure, the other at 500 psia. The functional 
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relationship for the head-loss ratio, equation 21, indicated 
that H* was dependent on void fraction, a, and pressure 
(for a given pump). In this equation, system pressure 
affects the two-phase flow function, a, through its effect 
on density and the slip-velocity ratio, s. Both of these 
effects are overshadowed by changes in void fraction, a. 

OTHER RESULTS 

The head-loss ratio was calculated using expressions 
other than equation 10 for the theoretical two-phase head 
coefficient. 

The experimental data for the fluid properties at 
pump-outlet were used to evaluate equation 9 and determine 
the head-loss ratio. Results are shown in Table El and 
Figure El. These results are not significantly different 
from the results presented in Figure 12, indicating that 
the effect of the outlet tangential velocity, Cq 2 > is 
minimal, as assumed. Equation 10 can be applied to any 
pump, without prior knowledge of the relationships between 
fluid properties at inlet and exit; therefore, it is much 
more useful as a means of predicting two-phase characteristics. 

Results obtained by assuming zero outlet-swirl are 
shown in Figure E2, and also show good correlation between 
H* and a, but indicate much higher losses. The procedure is 
simplified by this assumption, and further experiments may 
prove it to be viable for certain pump designs, particularly 
those of low specific speed. 

43 




4 


:« 


I 






APPLICATION OF THE PROCEDURE 


The two-phase, third-quadrant head-versus-flow 
characteristics of a centrifugal pump can be predicted, for 
various values of void fraction and pressure, using the pro¬ 
cedure outlined below. The following information concerning 
the pump must be known. 

Flow area at the pump discharge flange, A^• 

Flow area of the outer radius of the impeller, A^. 

Flow area at the impeller eye, A 2 . 

The impeller radius ratio, r 2 /r^. 

Blade angle at impeller eye, defined by 

Figure 3. 

Sufficient single-phase data to construct a plot 
of head-coefficient versus flow-coefficient. 

The actual single-phase, and the theoretical single- 
and two-phase head-versus-flow characteristics are determined 
by constructing a graph as shown in Figure 13. The three 

head-coefficients, ’^ithsp' '^ithtp' then be 

found for particular values of flow-coefficient, and 

two-phase flow function, 

At the desired value of void fraction, and system 

pressure, the head-loss ratio, H*, is determined from 
Figure 12, and the two-phase flow function, f^^^, calculated. 
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FIGURE 13 - ACTUAL SINGLE-PHASE AND THEORETICAL SINGLE- 
AND TWO-PHASE PERFORMANCE 
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The predicted values of actual two-phase head coeffi¬ 
cient, calculated as a function of flow-coefficient, 

(J)^, from 

♦itp = nthtp - »*>*ithsp 

Experiments on pumps of different geometries will be 
required to determine the effect of geometrical factors such 
as specific speed, N , on the functional relationship of the 
head-loss ratio, H* (equation 21). Until such data becomes 
available, this procedure may be used to determine approximate 
two-phase characteristics. 
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APPENDIX A 


DETAILS OF THE DERIVATION OF THE THEORETICAL 
TWO-PHASE HEAD COEFFICIENT 


In two-phase separated flow, the theoretical two-phase 
head coefficient is given by 


'^ithtp ^^“^1^^9Ll‘^^1^0vl^"y2 ^ ^2^^eL2‘^^2^ev2^ 


(Al) 


In general, velocities of the two phases are given by 






Pj. Aj. ' p^(l-a)A^ 


m 


V 


X 


C = -s-- 

V p A pa A_ 

^v V V T 


Assuming C^^ = C^ 


niT(l-Xi) 

^GLl ^ ^L3 ^ Pli {l-a^)A2 


= C 


^T 


'0vl v3 P^j^ A^ 


By continuity 


m,|, = constant 


^mLl Pli ^1 

1 - X, 


C . p , a, A^ 
mvl vl 1 1 

X, 
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Then 


^rriil ^Ll 


0L1 


(1-x^) 


Plj^(1-Oi)A3 


and 


9L1 


C — 
mLl 


Similarly 


C = C —i- 

evl mvl 


Substituting into equation Al 


C T 1 A, C , A, 

Wi ' = (1 —X ) —S^ii. —i 4- X _i. 

^Ithtp 1^ A^ 1 A^ 


(1 ^2^2^0v2, 

I o ■•■ o J 


U, 


U, 


From the general expressions for velocity 
^niLl ~ (l-Uj^) A^" 

and 


^mvl *^1 ^1 


(A2) 
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Substituting these expressions, the first two terms 


on the right hand side of equation A2 become 

2 • 2 • 

(1-x^) m^ x^ m^ A^ 

TT /T \ TV p a T UtA_ A- 

UiPli (l-cti) Ai 3 ^vl 1113 

Rearranging and multiplying by 


(1-a^) 


‘’lI * “l ‘’vl 


^Ll ^ “l 


where 


“'’v = '’tp 


They become 


“t 


2 2 
(1-x,)^ x; 

[t-f (l-aj^) p^j_) ] 


^1 '^tpl ^1 ^3 ^Ll^^ ‘^1^ ^vl “l' 


Substituting the definition of the two-phase flow coefficient 


m„ 


<P 


'^p'' Ptpi “i 


the two terms become 


a 7 ^tpi 


2 ^^"^1^ ^1 ^vl ^ ^1 ^L1 . _2 


Pvl ^‘l 


+ x^] 
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by definition, the two-phase flow function, a, is 


a = 


r_2L.] 

Pl 


and the quality, x, becomes 


X 


as 

1+as 


where 



substituting these expressions, the first two terms become 


A, 2a s^ 

^3 ^ ^ ^1®1 l+2a^s^+(a^s^) 


I (1-^ i) I 


performing the indicated operations and rearranging yields 


^1 . 

A7 ^tpl 


(1+a^) (l+a^^s^) 
(1+a^s^)^ 


Substituting this expression, equation A2 becomes 






Ithtp A^ tpl tpl 


, (1~X2)^2^0L2 , ^2^2^0v2t 
~ I o ■’* -) J 


U, 


u 
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where 


f 


tpl 


( 1 +a^) ( 1 +a^s^) 
(1+a^s^)^ 


From the outlet velocity triangle. Figure 4 , the last two 
terms in equation A2 reduce to 


[1 


^tp 2 *tp 2 . 2 

tan 62 


using the same procedure applied to the first two terms. 
Equation A2 becomes 


’^Ithtp A^ ‘^tpl ^tpl 


2 2 

-{^) [1 


^tp2 '*’tp2, 

tan B 2 


The two-phase flow coefficients, inlet and outlet 

can be related by the continuity equation. 


'^tp2 


'^tpl 


^tpl 

Ptp2 



Substituting this expression and rearranging yields 


ntKtp ' ^ ^PX 


+ fi *^tpl ^ ^tp2 j 

^tp2 ''2 ^2 


(A3) 
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APPENDIX B 


DETAILS OF THE CALCULATION OF HEAD-LOSS RATIO, H* 

The head-loss ratio, H*, for the various values of void 
fraction, a, was calculated using the relationships previously 
presented. The data from the experiments conducted by Olson 
are included in Appendix C. Values for fluid densities 
were interpolated from reference 9. 

Table B1 lists the single-phase characteristics used 
to determine the single-phase head coefficient, a 

function of flow coefficient, <}>^. Tables B2 through B5 
present the calculation of the head-loss ratio, H*. 

The test number column refers to the experimental data. 
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TABLE B1 


CALCULATION OF THE SINGLE-PHASE HEAD COEFFICIENT, ij;' , 



AS A FUNCTION OF FLOW 

COEFFICIENT 

* / 0 


SINGLE¬ 

PHASE 

TEST 

NUMBER 

gAH 

2 2 
(ft /sec) 

^1 

(ft/sec) 

^ml 

(ft/sec) 

Isp p 
(gAH/U^) 

^1 

78 

1330 

13.2 

3.188 

7.64 

0.24 

82 

1459 

14.2 

3.275 

7.22 

0.23 

86 

164.22 

13.4 

0.75 

0.91 

0.06 

87 

570 

13.5 

1.962 

3.15 

0.15 

88 

763 

13.5 

2.308 

4.22 

0.17 

89 

1346 

13.6 

3.188 

7.40 

0.23 

92 

14 

13.6 

0 

0.4 

0 

83 

4176 

20.8 

3.635 

9.63 

0.17 

75 

605 

27.1 

1.255 

0.82 

0.05 

76 

960 

27.2 

2.121 

1.30 

0.08 

77 

1513 

27.1 

3.231 

2.06 

0.12 

91 

299.5 

27.04 

0 

0.41 

0 

95 

2380 

53.78 

2.568 

0.82 

0.05 

98 

1204 

54.15 

0 

0.41 

0 

150 

4405 

54.25 

4.95 

1.50 

0.09 

106 

3130 

74.6 

1.414 

0.56 

0.02 

107 

3916 

74.7 

2.640 

0.70 

0.04 

108 

4634 

74.6 

3.534 

0.83 

0.05 

149 

7065 

74.5 

5.756 

1.27 

0.08 

147 

7467 

90.9 

5.02 

0.90 

0.055 

148 

9792 

91.2 

6.492 

1.18 

0.07 

152 

9982 

92.0 

6.65 

1.18 

0.07 
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TABLE B2 


CALCULATION OF TWO-PHASE FLOW FUNCTION, f , 

_ 


AND TWO-PHASE DENSITY, , AT PUMP INLET 


TWO- 

VAPOR 

LIQUID 

VOID 




PHASE 

TEST 

NUMBER 

DENSITY 

Pvl 

DENSITY 

3 

FRACTION 

Ui 

^tpl 

^tpl 

m 

149 

0.4377 

54.38 

0.548 


1.0299 

24.82 

150 

0.4443 

54.33 

0.461 


1.0218 

29.49 

151 

0.4508 

54.28 

0.464 


1.0151 

29.30 

152 

0.4354 

54.40 

0.466 


1.0218 

29.25 

153 

0.4354 

54.40 

0** 


1.0 

54.40 

154 

0.4420 

54.35 

0** 


1.0 

54.35 

155 

0.4309 

54.43 

0** 


1.0 

54.43 

156 

ic 

0.4288 

54.45 

0.241 


1.0150 

41.43 

157 

1.1002 

50.53 

0.116 


1.0086 

44.80 

158 

1.1066 

50.51 

0.245 


1.0069 

38.41 

159 

1.1042 

50.52 

0.056 


1.0039 

47.75 

160 

1.0932 

50.57 

0.512 


1.0671 

25.24 

161 

1.1156 

50.46 

0.537 


1.0232 

23.96 

162 

1.1156 

50.46 

0.631 


1.0327 

19.32 

163 

1.1453 

50.33 

0.401 


1.0143 

30.61 

164 

1.1202 

50.44 

0.647 


1.0348 

18.53 

*Tests 

149 through 

156 were 

conducted 

. at 

nominal 

system 

pressure of 200 psia, from : 

Figure 5, 

slip velocity ratio. 

s = 2.8 

. Tests 157 

through 

164 were 

conducted at 

nominal 

system 

pressure of 

500 psia 

, s = 2. 





* * 


Negative void fractions in data were taken to be zero 


54 






TABLE B3 


TWO- 

PHASE 

TEST 

NUMBER 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 


CALCULATION OF TWO-PHASE FLOW FUNCTION, f , 

- - 

AND TWO-PHASE DENSITY, AT PUMP OUTLET 


VAPOR LIQUID VOID 


DENSITY 

^2 

DENSITY 

Pl2 

FRACTION 

°‘2 

^tp2 

^tp2 

m 

0.4261 

54.47 

0.840 

1.1070 

9.07 

0.4113 

54.60 

0.867 

1.1234 

7.62 

0.4197 

54.52 

0.781 

1.0767 

12.27 

0.4197 

54.52 

0.558 

1.0298 

24.33 

0.4009 

54.68 

0.008 

1.0002 

54.25 

0.3967 

54.71 

0.007 

1.0002 

54.33 

0.3946 

54.73 

0.730 

1.0568 

15.07 

0.3864 

54.80 

0.789 

1.0742 

11.87 

1.0514 

50.76 

0.702 

1.0405 

15.86 

1.0646 

50.70 

0.756 

1.0509 

13.18 

1.0624 

50.71 

0.428 

1.0147 

29.46 

1.0646 

50.70 

0.696 

1.0400 

16.15 

1.0844 

50.61 

0.912 

1.1065 

5.44 

1.0844 

50.61 

0.988 

1.0860 

1.68 

1.1173 

50.46 

0.744 

1.0505 

13.75 

1.0998 

50.54 

0.834 

1.0736 

9.31 


(1+a)(1+as^) 

(1+as)^ 



(1-a) + ap^ 
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TABLE B4 


CALCULATION OF 

TWO-PHASE 

FLOW COEFFICIENT, 



AND TWO- 

PHASE HEAD 

COEFFICIENT, 

ntp 


TWO- 

PHASE 

TEST 

NUMBER 

FLOW 

RATE, Q 

(gal/min) 

ROTOR 
SPEED, N 

(rev/min) 

'^tpl 

(0.951 Q/N) 

HEAD 

AH 

o 

(ft) 

^itp 

(28,200 

2 ’ 
N^ 

149 

• 

in 

00 

1 

0 

— 

40.2 

— 

150 

-51.5 

-1608 

0.030 

84.3 

0.919 

151 

-30.3 

-1600 

0.018 

79.6 

0.876 

152 

-17.3 

-1607 

0.010 

74.8 

0.816 

153 

-16.2 

-1604 

0.009 

42.8 

0.469 

154 

-16.4 

-1602 

0.010 

43.8 

0.481 

155 

-39.8 

-1600 

0.024 

51.8 

0.570 

156 

-73.4 

-1603 

0.044 

75.1 

0.824 

157 

-80.0 

-1607 

0.047 

74.3 

0.811 

158 

-73.8 

-1603 

0.044 

77.3 

0.848 

159 

-42.6 

-1600 

0.025 

57.1 

0.629 

160 

-61.4 

-1594 

0.037 

90.5 

1.004 

161 

-36.2 

-1597 

0.022 

86.5 

0.956 

162 

-25.3 

-1600 

0.015 

107.5 

1.184 

163 

-26.5 

-1603 

0.016 

71.8 

0.788 

164 

-112.9 

0 

— 

76.4 

— 
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TABLE B5 


CALCULATION OF HEAD-LOSS RATIO, H* 


TWO- 

PHASE 

TEST 

NUMBER 

'^Ithtp 

*itp 

'^Ithsp 

■ Isp 

H* 

149 

— 

— 

— 

— 

— 

150 

0.095 

0.919 

0.089 

0.545 

1.81 

151 

-0.022 

0.876 

-0.025 

0.454 

1.88 

152 

-0.094 

0.816 

-0.096 

0.419 

1.77 

153 

-0.102 

0.469 

-0.102 

0.417 

1.10 

154 

-0.101 

0.481 

-0.101 

0.417 

1.12 

155 

0.027 

0.570 

0.027 

0.491 

1.17 

156 

0.215 

C.824 

0.209 

0.683 

1.28 

157 

0.247 

0.811 

0.244 

0.732 

1.16 

158 

0.214 

0.848 

0.211 

0.686 

1.33 

159 

0.043 

0.629 

0.042 

0.502 

1.27 

160 

0.168 

1.004 

0.146 

0.605 

1.82 

161 

0.012 

0.956 

0.008 

0.476 

2.02 

162 

-0.048 

1.184 

-0.052 

0.438 

2.51 

163 

164 

-0.043 

0.788 

-0.046 

0.442 

1.70 


nth = 

-0.1896 

+ 9.155 

^tpl "^l 



nsp = 

101.28 (p 

1.8756 

1 

0.4 



H* = ^ithtp ~ "^itp 


INLET 

VOID 

FRACTION 


0.461 

0.464 

0.466 

0 

0 

0 

0.241 

0.116 

0.245 

0.056 

0.512 

0.537 

0.631 

0.401 
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APPENDIX C 


EXPERIMENTAL DATA 

Tables Cl and C2 are reproduced from reference 7. 

Explanatory notes for single-phase tests: 

a. Based on turbine flowmeter measurements. 

b. Positive pump AP represents higher pressure 
at discharge nozzle. 

c. Based on fluid density measurement at pump 
inlet. 

d. Fluid density at pump inlet (suction for 
normal operation, discharge for reverse flow). 

e. Inlet and outlet refer to locations with 
respect to direction of flow. 

f. Flow rate obtained from venturi flowmeter 
in test loop. 

Explanatory notes for two-phase tests: 

a. Inlet flow based on loop mass flow rate and 
test section inlet (suction or discharge) 
fluid density. 

b. Positive AP indicates higher pressure at 
pump discharge. 

c. Based on inlet conditions. 

d. Inlet and outlet refer to the direction of 
flow. 
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e. Calculated using measured density and 
assumed saturation conditions. 

f. Measured directly. 
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RESULTS OBTAINED FROM TEST SECTION MEASUREMENTS DURING SINGLE-PHASE PLOW PUMP TESTS 
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APPENDIX D 


DESCRIPTION OF EXPERIMENTS 

The semi-scale pump used in the experiments had the 
following characteristics: 

Rated flow - 180 gpm 
Rated head - 192 ft. 

Specific speed - 926 
Rated torque - 417.6 in-lb 
Details of the pump are shown in Figures 9 and 10. 

The test-loop used for single-phase tests is shown in 
Figure Dl. The two-phase test-loop is shown in Figure D2. 
Superheated steam and/or saturated water were supplied by 
an oil-fired boiler. 

Flow conditions imposed for third quadrant tests are 
as follows. 

Single-phase tests 

Impeller speed - up to minus 2700 rpm 
Flow rate - 0 to minus 207 gpm 
Pressure at pump inlet - 60.7 to 216.6 psig 
Pump head - 5.1 to 310.0 ft of fluid 
Two-phase tests 

Impeller speed - 0 and minus 1600 rpm 

Flow rate - minus 16.2 to minus 112.9 gpm 

System pressure - 200 and 500 psia 

Pump head - 40.2 to 107.5 ft of fluid 
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FIGURE D1 - SINGLE-PHASE TEST LOOP 
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(from reference 7) 
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FIGURE D2 - TWO-PHASE TEST LOOP (from reference 7) 
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APPENDIX E 


OTHER RESULTS 

This section summarizes, for purposes of comparison, 
the correlation of H* vs a obtained using alternate methods 
to compute the theoretical characteristics. 

The experimental data were used to evaluate equation 9, 
thereby making it semi-empirical rather than theoretical. 

The head coefficients obtained in this manner were used to 
calculate the head-loss ratio. Results are shown in Table El 
and Figure El. 


ZERO SWIRL AT OUTLET 

The theoretical characteristics can be derived assuming 
that outlet swirl is negligible. In this case the expression 
for the theoretical head coefficient, is 


^1 

’^ith A^ ^tpl '*’1 


Results are shown in Table E2, and the correlation of 
H* versus is shown in Figure E2. 
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TABLE El 


CALCULATION OF THE HEAD-LOSS RATIO, H*, USING 
EXPERIMENTAL DATA TO EVALUATE EQUATION 9 


TWO- 

PHASE 

TEST 

NUMBER 

^'^ithtp^ se 

H* 

°‘l 

150 

0.199 

1.581 

0.461 

151 

0.006 

1.816 

0.464 

152 

-0.092 

1.763 

0.466 

153 

-0.102 

1.100 

0 

154 

-0.101 

1.124 

0 

155 

0.096 

1.023 

0 

156 

0.336 

1.030 

0.241 

157 

0.341 

0.964 

0.116 

158 

0.306 

1.142 

0.245 

159 

0.059 

1.238 

0.056 

160 

0.189 

1.776 

0.512 

161 

0.097 

1.835 

0.537 

162 

0.128 

2.156 

0.631 

163 

-0.022 

1.660 

0.401 

'^Ithtp 

= -0.1896 + (8.1342 

f^ , + 1.021 
tpl 

^tpl ^ 


values for ftp2' '^tpl' ‘^tp2 obtained from 

Appendix B. 


t 


where subscript se 


semi-empirical 


69 

























































































FIGURE E2 - HEAD-LOSS RATIO VERSUS INLET VOID FRACTION FOR 
ZERO-SWIRL AT OUTLET, - 0* 
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TABLE E2 



CALCULATION OF THE 

HEAD-LOSS RATIO, 

H*, 



ASSUMING ZERO SWIRL 

AT ROTOR OUTLET, 

^0 O . 

= 0 

TWO- 

PHASE 

TEST 

NUMBER 

’^Ispth 

'^itpth 


H* 

150 

0.248 

0.253 


2.24 

151 

0.146 

0.148 


2.36 

152 

0.083 

0.085 


2.18 

153 

0.078 

0.078 


1.15 

154 

0.079 

0.079 


1.19 

155 

0.193 

0.193 


1.27 

156 

0.354 

0.359 


1.41 

157 

0.385 

0.388 


1.22 

158 

0.356 

0.358 


1.48 

159 

0.206 

0.207 


1.43 

160 

0.298 

0.318 


2.23 

161 

0.175 

0.179 


2.58 

162 

0.122 

0.126 


3.35 

163 

0.128 

0.130 


2.10 
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